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ABSTRACT 

We present XMM-Newton RGS and EPIC data of the putative cooling flow cluster Abell 2597. 
Velocities of the low-ionization emission lines in the spectrum are blue-shifted with respect to the 
high- ionization lines by 1 320i®]o kms -1 , which is consistent with the difference in the two peaks 
of the galaxy velocity distribution and may be the signature of bulk turbulence, infall, rotation 
or damped oscillation in the cluster. A hierarchical velocity structure such as this could be the 
direct result of galaxy mergers in the cluster core, or the injection of power into the cluster gas 
from a central engine. The uniform X-ray morphology of the cluster, the absence of fine scale 
temperature structure and the random distribution of the the galaxy positions, independent of 
velocity, suggests that our line of sight is close to the direction of motion. These results have 
strong implications for cooling flow models of the cluster Abell 2597. They give impetus to those 
models which account for the observed temperature structure of some clusters using mergers 
instead of cooling flows. 

Subject headings: galaxies: clusters: individual: A2597 - galaxies: distances and redshifts - X-rays: 
galaxies 


1. Introduction 

It has long been theorized {Burns et al. 2002) 
that the formation of clusters by hierarchical 
mergers would produce velocity substructure in 
the X-ray emitting gas in clusters of galaxies. 
However the observational consequences of such 
substructure is not clear and several detailed sim- 
ulations differ in their predicted amplitude and 
detailed structure (c.f. Roettiger k Mushotzky, 
1998; Norman k Bryan 1999; Nagai k Kravtsov 
2003; Onuora, Kay and Thomas 2003). The pre- 
dicted structures range from generalized turbu- 
lence, the existence of large velocity structures, 
the signature of the infall of an individual group, 
rotation and infall along the filaments which meet 
at the cluster. Despite this uncertainty, the gen- 
eral existence of velocity structure is a strong pre- 
diction of hierarchical merging. The other possi- 
ble source of velocity structure is the action of an 
AGN in the center of the galaxy which can put me- 
chanical energy into the gas via the production of 
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jets, radio “bubbles” or other unspecified means, 
producing bulk motion and turbulence (e.g., Chu- 
razov et al. 2000, Reynolds, Heinz and Begelman 
2002, Fabian et al. 2003). 

Thus, while long predicted, the only evidence to 
date for the existence of such motion is the sup- 
pression of resonance scattering in several clus- 
ters and galaxies (Xu et al. 2002; Gastaldello 
k Molendi 2004; Churazov et al. 2004) and the 
shapes of head-tail cluster radio sources (Burns 
et al. 2002). The relative level of resonance scat- 
tering suppression implies fairly strong turbulence 
or velocity shear with a Mach speed > 1/3 of the 
sound speed in the X-ray emitting gas associated 
with the radio sources in M87 and the Perseus 
cluster. However in the only data so far pub- 
lished of a radio quiet giant elliptical galaxy, NGC 
4636, there is a strong upper limit on turbulence 
or velocity shear (Xu et al. 2002). There is, at 
present, no strong evidence for velocity structure 
in the X-ray emitting gas in other clusters, groups 
or galaxies (but see Bregman and Dupke 2002 
for possible large scale velocity structure in the 
Perseus and Centaurus clusters). 
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In this paper we report the detection with the 
XMM-Newton RGS gratings of velocity substruc- 
ture in the core of of the rich cluster Abell 2597, 
one of the “massive cooling flow clusters” (Craw- 
ford et al. 1989). This velocity structure is consis- 
tent with the observed redshift distribution of the 
galaxies. Rather than a cooling flow, we believe 
that this data supports the merger model (Motl 
et al. 2004) for the origin of the cool gas in this 
cluster. 

2. Observations 

XMM-Newton (Jansen et al. 2001) observed 
Abell 2597 between 03:29 UT and 09:20 UT on 
Nov 30, 2000 2 . EPIC, RGS and OM cameras were 
in operation; see Striider et al. (2001), Turner 
et al. (2001), den Herder et al. (2001) and Mason 
et al. (2001) for instrument details. Raw Obser- 
vation Data Files files were reprocessed using the 
pipeline software distributed within SAS v5.4 and 
spacecraft calibrations which were up-to-date as 
of Dec 11, 2003. Since there was no significant 
proton-induced background flare activity during 
the exposures, additional time-screening was not 
required. 

The EPIC MOS and EPIC pn cameras were op- 
erated in full-frame mode behind medium filters. 
Event files were screened for known and new bad 
pixels, and events of complexity greater than dou- 
bles were rejected. Events with corrected pulse 
heights below 300 eV were screened out to avoid 
calibration uncertainties. 

RGS first order spectra were extracted from 
a source position of RA=351. 33229 and Dec=- 
12.12389 using channel binning of 0.045A. The 
cross-dispersion extraction region was defined to 
include 98 percent of the flux within the RGS 
cross-dispersion Point Spread Function (PSF). 
The shape of the cross-dispersion PSF is variable 
over the dispersion axis, with the extraction region 
ranging in width between 2.2-3.8 arcmin. With a 
5 arcmin field of view in the cross-dispersion di- 
rection, this allowed sky regions to be extracted, 
albeit with reduced signal-to-noise compared to 
point source data and with some small contami- 
nation from the faint, outer regions of the cluster. 
Due to the small number of events, second order 
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spectra are not analyzed. 

3. EPIC results 

We performed imaging spectral analysis of 
Abell 2597 in the central 6 arcmin, where the sur- 
face brightness exceeds the internal background. 
Extracting events from concentric annuli about 
the cluster center and fitting collisional plasma 
models (Brickhouse et al. 2000), we measure an 
average overall temperature of 3.5 keV consistent 
with the ASCA average temperature and evidence 
for a temperature drop in the central 50 arcsec, 
consistent with the Chandra results of McNamara 
et al. (2001). While our temperatures are statisti- 
cally more precise than those of McNamara et al. , 
the poorer angular resolution does not allow full 
resolution of the central temperature profile. The 
X-ray color map shows little evidence for spatial 
temperature structure other than the overall re- 
duction in temperature in the central 45 arcsec. 
The source was found to be isothermal at radii 
> 1 arcmin from the center with no abundance 
gradients. In Fig. 1 we plot the temperature 
map, with intensity contours overlayed, and the 
best-fit temperatures and Fe abundances versus 
angular distance from the cluster core. 

4. RGS results 

First order RGS spectra of Abell 2597 are plot- 
ted in Fig. 2, revealing emission lines of O vm 
Lya, L-shell Fe transitions and, marginally, Si xiv 
Lya. 

4.1. Emission line velocities 

Because the target source is extended, there is a 
range of incident angles on the gratings and emis- 
sion lines are consequently broadened. In order 
to fit the line profiles adequately we must simu- 
late the instrumental broadening. The simplest 
method is to convolve a model spectrum with the 
integrated spatial profile of the source collapsed in 
the cross-dispersion direction of the RGS cameras, 
producing a convolution function for the model. 
This is achieved with the model rgsxsrc , added to 
the XSPEC package for version 11.2. In this case we 
have used a 0.35-2.5 keV image, extracted from the 
EPIC MOS1 events list to generate the function. 
Although fast computationally, the disadvantage 
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Fig. 1. — Left: X-ray hardness ratio image of Abell 2597, extracted from the EPIC MOS event data. White 
to red indicates increasing temperature. The bands employed in the calculation were 0.3-2 keV and 3-8 ke\ 
with images smoothed by a factor ??? using Gaussian convolution. Counts of 6-8 per pixel correspond to 
blackbody temperatures of 3. 6-4. 6 keV, 8-10 counts to 2. 9-3. 6 keV and 10-12 counts to 2. 5-2. 9 keV. X-ray 
intensity contours from an unsmoothed 0.3-2 keV image combined from MOS1 and MOS2 are superposed 
over the image. Right: Collisional plasma fits to EPIC spectra extracted from concentric annuli around the 
cluster core. We plot plasma temperature and Fe abundance as functions of angular distance from the core. 
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Wavelength (A) 

Fig. 2.— XMM-Newton RGS1 (blue) and RGS2 (red) spectra of Abell 2597. The black lines correspond to 
the best-fit two-temperature collisional plasma model, described in Sec. 4. Lines from 0 vm, Fe xxm (or 
Fe xxn), Fe xxiv, and possibly Si xiv are detected. The bottom panel plots the fit residuals. 


of this method is that the emitted spectrum is as- 
sumed to be isothermal across the source. 

To measure the redshift of the O vm Lya line 
we use a Gaussian model, with a powerlaw to rep- 
resent the continuum across the wavelength range 
19-22A. The best fit model, convolved with the 
spatial profile above, is plotted in Fig. 3. We take 
the rest-energy of the line as 653.6 eV (Johnson 
& Soff 1985) and assume that the local intrinsic 
width of the line is negligible. The redshift of this 
feature is consequently z = 0.08lj^oo4 (x 2 = 51 
for 52 d.o.f.), where the errors are 90 percent confi- 
dence limits. This is less than the commonly cited 
mean optical redshift of the cluster z = 0.0852 
(Struble and Rood 1999) with greater than 99 per- 
cent confidence. 

An identical model, ray-traced through a 
simulation of the XMM-Newton telescope us- 
ing the Monte Carlo code discussed in Peter- 
son et al. (2001), yields a consistent fit of z — 
0.078 ± 0.001. While the result is not particu- 
larly sensitive to the spatial distribution of the 
source, based on the broad band EPIC MOS im- 
age, we assumed the O vm line is distributed as a 
/3-function, where the spatial profile has an index 


of 0.62 and core radius of 27 arcmin. 

A convolved Gaussian plus powerlaw model fit 
to the suspected Lya Si XIV line over the AA5- 
8A range provides z = 0.078lo oo9 (X 2 = 52 for 
52 d.o.f.), taking the rest wavelength as 2.005 keV 
(Johnson & Soff 1985). This is consistent with the 
O vni redshift and, again, inconsistent with the 
Struble and Rood redshift with 90 percent confi- 
dence. 

To fit the Fe L shell transitions we adopt a 
neutrally-absorbed thermal plasma model ( vapec ; 
Brickhouse et al. 2000) over AA8-15A, con- 
volved with the EPIC MOS1 spatial profile. 
Photo-electric cross-sections were provided by 
Balucinska-Church & McCammon (1992). Best fit 
parameters are z = 0.086lo_oo5> n H = ^^xlO 20 
cm -2 and kT — 3.ltJ;g keV, with an Fe abun- 
dance of 0.44t°;i9 solar (x 2 = 189 for 179 d.o.f.). 
The velocity of the Fe emission is consistent with 
both the Struble and Rood optical redshift of the 
cluster and the O VIII line within the 90 percent 
confidence limits. However there is one strong 
line at A12.7A that is not fit well by the model. 
This is partly because it is both blue-shifted rela- 
tive to the best-fit velocity and narrower than the 
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Fig. 3. — RGS1 (blue) and RGS2 (red) spectra of the 0 vm line from Abell 2597. The black lines correspond 
to a powerlaw + Gaussian fit, convolved with the 0.35-2.5 keV EPIC MOSl spatial profile. The redshift of 
the line corresponds to z = 0.081, whereas the mean optical redshift is represented by the vertical dotted 
line. The lower panel plots fit residuals. 




EPIC MOSl convolution function. An additional 
line with negligible intrinsic width was added to 
the model and fit to this feature. It is likely this 
line is either Fe xxn A11.770A at z — 0.0801q^q° 4 
or Fe xxm All. 736.4 at z = 0.082lg;ggJ (Brown 
et al. 2002). 

The RGS2 chip gap between AA13.1-13.2A may 
overlap another strong line. The current model is 
fitting solar abundance Ne x Lycc to the red wing 
of the feature, but this line equally may be Fe XXI 
A 12. 284 A, which would reduce the temperature of 
the model considerably. By removing Ne from the 
model completely, best fit thermal plasma param- 
eters are z = 0.085±g:gg|, n H = 9tf xlO 20 cm' 2 
and kT — 2.0to'es keV, with an Fe abundance of 
0.22t°;Jo solar (x 2 = 187 for 179 d.o.f.). 

The temperature is still too high to produce 
significant amounts of Fe xxi (Arnaud & Rothen- 
flug 1985), but this is not especially surprising. In 
reality, the source has a range of plasma temper- 
atures (Sec. 4.2), while the fit above converges 
on a mid-range value. The quoted tempera- 
ture is also artificially high because the narrow 
Fe xxn/Fe xxm line is not being fit correctly by 
the instrumentally-broadened model. What we 
can say with some certainty is that in collisional 
conditions, the absence of ionization states with 
energies equal or less than Fe xx indicates plasma 
temperatures > 1 keV, while if Fe xxn is the cut- 
off state, the plasma temperature must be > 1.5 
keV (Arnaud & Rothenflug 1985). 

In summary, the hot Fe xxiv lines have veloci- 
ties consistent with the Struble & Rood (1999) op- 
tical redshift of the cluster. The unblended 0 vm 
and Si xiv Lycc lines are blueshifted relative to 
this velocity with > 90 percent confidence. The 
velocity of the cooler Fe transition, either Fe XXII 
or Fe xxm, is statistically more consistent with 
the O VIII and Si xiv and is narrower than the 
Fe xxiv lines, suggesting an origin within a local- 
ized region of the cluster. 

4.2. A two-zone model 

Fe transitions of different widths suggest that 
these lines are spatially distributed differently in 
the source. The Fe line widths are as expected 
phenomenologically for a cooling flow - the broad 
lines arise in a large, hot outer volume while the 
narrow lines originate in the cooler, inner flow. 


This model explain explains elegantly why the 0 
and Si lines are broad, yet display the peculiar 
velocity of the cool core. Whereas the Fe XXII 
and Fe xxm ionization states emit over relatively 
narrow distributions in temperature, H-like 0 and 
Si have distributions considerably broader, with 
ranges encompassing Fe xxi-xxiv. Consequently, 
the 0 and Si lines should arise in both the inner 
and outer regions of the source. 

A first order approximation to this picture is a 
two-zone model. The outer zone is an extended 
thermal plasma which we model using the EPIC 
MOSl convolution function to broaden the spec- 
tra. The inner zone is a thermal plasma, mod- 
eled as a point source. Both zones are absorbed 
by the same galactic neutral hydrogen column in 
the direction of the source, tih = 2.49 xlO 20 cm -2 
(Dickey & Lockman 1990). 

Plasma parameters for the outer zone were fixed 
at values determined from the EPIC data (Sec. 
3). We take kT 0 = 3.5 keV and element abun- 
dances [Fe/H] = —0.52 and [O/H] = —0.32. The 
suffix “o” corresponds to the outer region of the 
model. Since the source was found to be isother- 
mal at radii > 1 arcmin from the center with no 
abundance gradients, our choice of an isothermal 
convolution template is vindicated. Free floating 
parameters for this component are the normaliza- 
tion and the Si abundance which, being poorly 
constrained by the EPIC data, is tied to the Si 
abundance of the point source component in order 
to reduce redundancies. The temperature, kT \ , Fe, 
O and Si abundances, and normalization of the in- 
ner zone component are also free floating param- 
eters. The velocities of the two zones are initially 
coupled, yielding a best-fit redshift for the clus- 
ter of z — 0.083+Q OO4 (x 2 = 185 for 179 d.o.f.). 
This value is not consistent with the Struble and 
Rood (1999) optical redshift within the 90 percent 
confidence limits. Furthermore, the formal chi 2 is 
driven mostly by optically thin continuum which 
is insensitive to the redshift. This current fit rep- 
resents the velocities of the O vm and Fe xxill 
lines poorly. 

The fit to these lines improves if we decouple 
the redshift of the two zones. First, by fixing 
the redshift of the outer zone to the Struble and 
Rood value, z 0 = 0.0852. The best-fit inner zone 
redshift in this case is z, = O.O8O8I0 0022 (x 2 = 
184 for 179 d.o.f.). Best fit parameters for this 
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model are provided in Table 1, and the best fit 
model spectrum is plotted over the data in Fig 2. 
z\ = 0 . 0808^0 0022 * s l ess than the mean optical 
redshift with 8cr confidence. Conversely if we fix 
zi at 0.0808, z 0 = 0.0901°;°°® ( X 2 = 183 for 179 
d.o.f.) 

For comparison, a two zone model applied to 
the central 1 arcmin of the EPIC data gives a best 
fit temperature of kT; = 1.7i{j;J keV, perfectly 
consistent with the RGS model. The relative emis- 
sion measure ratio of the hotter to the cooler com- 
ponent dropping from 6:1 at angular distances R 
— 30-60 arcsec to 1.3:1 inside 15 arcsec. There 
is little evidence for cool, kT < 3 keV, gas out- 
side 90 arcsec in the EPIC spectra. In these two 
component fits there is no change in the tempera- 
ture of the cooler component with radius. Fitting 
a cooling flow model (model mkcflow in xspec) 
to the same regions, we find that the minimum 
temperature in the cooling flow is kT ~ 1.6 keV, 
as is often found in fits of cooling flow models to 
XMM-Newton and Chandra data. 

5. Discussion 

We have found evidence that the low temper- 
ature X-ray emission lines of O vm, Si xiv and 
Fe xxill from Abell 2597 have a different redshift 
to the high temperature lines of Fe xxrv. In this 
section we exclude the possibility of systematic 
bias to the redshift measurements, compare our 
results to data at other wavelengths, and discuss 
these findings in the context of cooling flows, group 
mergers and the influence of a central, active nu- 
cleus. 

5.1. Pointing uncertainties 

In this section we briefly discuss whether the 
measured RGS velocities can be affected by tele- 
scope pointing bias. There is a small discrep- 
ancy between the location of the cD galaxy in 
the Digitized Sky Survey (spatially consistent with 
the brightest knot in the high resolution Chandra 
ACIS image; Fabian et al. 2001) and the position 
of the source in the OM images. This is most 
likely a symptom of the telescope pointing accu- 
racy, which is limited to ~ 5 arcsec. With a po- 
sition angle of 247 degrees, the disparity is ~ 0.4 
arcsec in the dispersion direction of the RGS de- 
tector. To determine what effect an off-axis point- 


ing on these scales will have on the point-source 
wavelength scale, we consider the derivative of the 
dispersion equation: 

AA = dsinaoAa (1) 

using a nominal incidence angle of ao = 1.5762 de- 
grees and a grating groove density of d~ l = 645.6 
lines mm -1 (den Herder 2001). After correction 
for optical geometry, the angular discrepancy is 
An = 0.45 arcsec. This yields a drift in the wave- 
length scale at Ovm Lya of AA = 1 mA. The 
redshift difference of our two model components 
at the wavelength of the Ovm line corresponds 
to A A ~ 0.1 A. For pointing inaccuracies to con- 
tribute significantly to our measurements, we re- 
quire a pointing offset of order 48 arcsec. There- 
fore the effect of pointing uncertainties are negli- 
gible. 

5.2. Cluster velocity distribution 

It is constructive to compare the two measured 
RGS velocities above with the redshift distribu- 
tion of Abell 2597 cluster members. The often 
cited redshift of the cluster, z = 0.0852, is taken 
from the catalogue of Struble and Rood (1999) 
and it originates from the optical radial velocities 
of only three cluster members, corrected to the Lo- 
cal Group standard of rest (Kowalski, Ulmer and 
Cruddace 1983). None of these sources are the cD 
galaxy. To similarly correct our own results to the 
Local Group, we must add 140 kms -1 (Courteau 
and van den Bergh 1999). This removes just 10 
percent of the discrepancy between the cool X- 
ray lines and the Struble & Rood (1999) redshift. 
The standard deviation between the three galac- 
tic velocities, z = 2x 10 -3 , is of the same or- 
der as the difference between our two RGS ve- 
locities, but also of the same order as the sys- 
tematic uncertainty introduced by using different 
cross-correlation templates in the optical redshift 
measurements. 

The recent availability of the the 2dF galaxy 
survey data has allowed a drastic improvement in 
the measurement of the redshift and distribution 
of galaxy velocities (De Propris et al. 2002). Us- 
ing the latest 2dF data release, we find 44 galaxies 
within 30 arcmin of the center of Abell 2597 and 
their redshift distribution is shown in Fig. 5. It is 
immediately clear that there axe 2 redshift peaks 
which correspond to the velocities seen in the RGS 
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Table 1 


Best fit parameters for the inner, point-like component of the two-zone emission model 
for Abell 2597. Outer zone parameters have been fixed using EPIC measurements and the 

REDSHIFT OF THE OUTER ZONE COMES FROM STRUBLE AND ROOD (1999). 
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Fig. 5. — Histogram of the galaxy redshifts available from the 2dF survey in the Abell 2597 field. Two peaks 
in the redshift distribution are found, z ~ 0.082, consistent with the cool component of the RGS model, and 
z ~ 0.084, consistent with the hot component. The dashed line corresponds to the redshift of the cD galaxy, 
z — 0.0822 (Owen, Ledlow and Keel 1995). 


8 



data. The lower redshift peak is the one published 
by De Propris et al. (2002) and the higher peak is 
the one cataloged by Struble and Rood (1999). 
The lower redshift system is also consistent with 
the velocity of Ovi (Oegerle et al. 2001), the ra- 
dial velocity of the cD galaxy itself of z = 0.0822 
(Owen, Ledlow and Keel 1995), the velocity of 
Hi seen in absorption against the radio source 
(O’Dea, Baum and Gallimore 1994), and the ve- 
locity of Lya seen in the archival STIS data. We 
make a cut in the galaxy redshift distribution at 
z = 0.083 and throw away sources with z < 0.06 
or z > 0.1. Plotting the positions on the sky of the 
two velocity groups in Fig. 6, we find no segrega- 
tion in right ascension or declination. This indi- 
cates that the velocity structure is approximately 
orthogonal to the plane of the sky. 

5.3. Interpretation 

We report, for the first time, significant veloc- 
ity differences between ion species in the hot X-ray 
emitting gas in a cluster. The velocity differences 
between the lines produced in gas with kT <2 keV 
and the Fe xxrv lines produced at higher effective 
temperatures is Az = (4.4tl'%) x 10 -3 or 1 320^210 
kms -1 , and is significant at > 99 percent confi- 
dence. This velocity separation is completely con- 
sistent with the observed bi-modal galaxy veloc- 
ity distribution (Fig. 5), and the lower- velocity, 
cooler gas lies at the redshift of the cD galaxy. 
This range in gas velocity is slightly supersonic 
since the sound speed is ~ 950 kms -1 for the 
kT = 3.5 keV gas which constitutes the bulk of 
the cluster emission. 

At shorter wavelengths, the Hi absorption spec- 
trum towards Abell 2597 shows complex velocity 
structure over the central 3 axcsecs (O’Dea, Baum 
and Gallimore 1994), similar to that detected in 
the Perseus cluster. Projected across the nucleus 
is a clump, infalling by ~ 300 kms -1 (as opposed 
to our outflow of 1320 kms -1 ) and a broad Hi 
absorption (~ 400 kms -1 ) relative to the sys- 
temic velocity of the cD galaxy at z = 0.082. The 
Ha emission (Heckman et al. 1989) has a velocity 
width of ~ 600 km s -1 , but is also at the systemic 
velocity. Thus while there is strong evidence of 
velocity turbulence and shear in Abell 2597 from 
radio and optical data, neither reveal evidence for 
redshifted gas. There is also no evidence in the 
Ovi or Lyman alpha emission for the higher red- 


shift component. We thus associate the lower red- 
shift gas with the velocity of the cD and its asso- 
ciated galaxies. 

5.3.1. Cluster merger 

The cool gas velocity relative to the hot gas 
is quite consistent with with the observed galaxy 
distribution and projected velocities predicted in 
many numerical simulations. As first reported in 
Roettiger, Burns and Loken (1996), the gas in the 
lower temperature merging groups remains bound 
to the group for a considerable length of time. 
Thus it is possible that a “cool” core is, in reality, a 
projection of a merging group as it falls into a rich 
cluster. If this is so, this gas would be expected to 
retain the velocity of its infall for a considerable 
period of time. This scenario is discussed in de- 
tail by Motl et al. (2004), who predict that many 
of the “cooling flows” observed are in reality pro- 
jections of infalling groups. To first order in this 
model, one predicts velocity shear across the cen- 
ter of the cluster due to projection effects, where 
the observed cool core is due to emission of an in- 
falling group, projected against the cluster center. 
This idea is also consistent with the relative lack of 
cooling gas seen in the XMM-Newton RGS spectra 
of many clusters (Peterson et al. 2003). However, 
as stressed by Motl et al. , while the predicted ve- 
locity structure is quite complex, material is al- 
ways streaming through the cluster core. Thus 
the Motl et al. results are not well described by 
our simple model of two components to the veloc- 
ity structure. However it is clear that the merging 
model for cooling flows naturally predicts the ve- 
locity structure seen in our analysis and allows a 
natural interpretation of the data. 

There are reasons to doubt this interpretation. 
The X-ray image and temperature map (Fig. 1) 
show little structure with the X-ray isophotes be- 
ing slightly elongated in the NW-SE direction, 
consistent with the galaxy distribution direction. 
We have quantified this statement by extract- 
ing spectra from regions of the X-ray color im- 
age that deviate by more than la from the mean 
color and with more than 2000 counts in the 
PN spectrum. We find no evidence for temper- 
ature structure beyond that given by statistical 
error, with each region having internal errors of 
approximately ±0.5 keV. Both the XMM-Newton 
and Chandra (Fabian et al. 2001) azimuthally av- 
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Fig. 6. — Positions of sources within the 2dF survey have been plotted over a digitized plate from the 
Palomar sky survey. Source with z < 0.06 or z > 0.1 have been omitted. Remaining galaxies have been 
segregated into two redshift bins, where the division is made at z = 0.083. The two populations do not have 
different distributions on the sky. The cD galaxy is listed in the 2dF survey and is centered in the image. 
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eraged temperature profiles show that there is a 
significant reduction in cluster temperature only 
in the central 1 arcmin. Thus if there is an in- 
falling group it has to be rather small, < 1 ar- 
cmin in projected size (~ 100 kpc), and well 
centered on the overall X-ray surface brightness 
isophotes. While this seems rather unlikely, the 
probability of detecting a purely radial inflow of 
a group along a filament is not well known. Nor- 
man and Bryan (1999) estimate using the Lacey 
et al. (1993) merger tree that an appropriate mi- 
nor merger occurs every 0.4 Gyr, and that the ve- 
locities decay on a similar time scale. Since cool 
cores in clusters are very common, this seems quite 
reasonable, but the relative solid angle of inflows 
along filaments has not been well determined. 

The observed galaxy distribution within Abell 2597 
(Trevese et al. 1992) is highly elongated in the 
plane of the sky (Fig 7), consistent with seeing a 
filament at a large inclination angle. We thus spec- 
ulate that the combination of the X-ray velocity 
structure, galaxy velocity distribution, lack of X- 
ray temperature structure and galaxy distribution 
can be explained by a fortunate set of circum- 
stances in which we are observing a merger of two 
objects at a very oblique angle to the plane of the 
sky, and in which there are no strong shocks in 
the merger. However a large compilation of sim- 
ulated clusters (http://sca.ncsa.uiuc.edu) shows 
that such occurrences may be rare. 

5.3.2. Central engine 

Alternatively the velocity of the cool central 
region could be due to the effects of an active 
galaxy. In most of the Chandra spectral images 
of X-ray emission influenced by the central radio 
source (e.g. Belsole et al. 2001; Fabian et al. 2002; 
Blanton, Sarazin and McNamara 2003), the re- 
gion around the radio emitting plasma is cooler 
than the surrounding gas. The central galaxy of 
Abell 2597 (Sarazin et al. 1995) is a luminous ra- 
dio source and there are indeed ’’cavities” in the 
X-ray emission similar to those seen in Perseus 
and Abell 2052 (McNamara et al. 2001). However, 
analysis of the Chandra hardness ratios shows that 
the emission immediately adjacent to each cav- 
ity is generally no harder or softer than its sur- 
roundings (McNamara et al. 2001) . In addition, 
in neither the Chandra nor XMM-Newton images 
is there evidence for a central point source indi- 


cating that the present day luminosity of a central 
AGN is rather low. As noted in McNamara et al. , 
if the X-ray cavities are due to buoyant bubbles, 
they must rise in the central regions on timescale 
of 10 7 ' 7 — 10 8 yrs, and thus have a velocity of 300- 
500 kms -1 . This is somewhat slower than our 
observed velocity. However, since we see the bub- 
bles projected in the plane of the sky, the true 
velocity may be larger. 

In the numerical simulations of Reynolds 
et al. (2002) and Bruggen (2003) one clearly sees 
velocity structure still present long after the radio 
source has turned off. In the Omma et al. (2003) 
simulations, velocities of up to 4000 kms -1 are 
obtained by X-ray emitting gas along the axes of 
jets, even at distances of ~ 50 kpc from the center. 
It is not clear in this scenario what the origin of 
the single velocity shift we observe is, since one 
expects more velocity broadening than velocity 
shifts. However as these authors point out, about 
17 Myr after the jet turns off, the vortex fills and 
overtakes the cavity, which is thus replaced by 
a region of outward-moving, cool, overdense gas, 
which may produce our observed outgoing velocity 
structure. 

An alternative source of high velocity gas is the 
copious star formation occurring in the central re- 
gions of Abell 2597 (Koekemoer et al. 2002), in 
which case we may be observing the equivalent of 
a ’’blowout” in a starburst galaxy. However, since 
the gas is expanding into a hot dense environment 
in this case, the boundary conditions are rather 
different and it is not clear what the velocity field 
should be. 

It is tempting to speculate, based on the Chan- 
dra and XMM-Newton results for the Perseus clus- 
ter (e.g. Fabian et al. 2003, Churazov et al. 2004; 
Gastaldello and Molendi 2004), and the Chan- 
dra images of Abell 2597 (McNamara et al. 2000), 
that our observed velocity structure is due to tur- 
bulence in the inter-galactic medium, related to 
the presence of the radio structure. However this 
would then not explain the coincidence with the 
velocity structure seen in the galaxies. 

5.4. Effects on cooling flow models 

The effect of velocity structure on models of 
cooling flows is severe. Abell 2597 had an inferred 
cooling rate of ~ 300-400 M 0 per year (Crawford 
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Fig. 7. — A contour map of the XMM-Newton EPIC M0S1 image overlayed on a density map of galaxies. 
Sources were identified from the Palomar survey by Trevese et al. (1992). We smoothed the galaxy map by 
a Gaussian function of FWHM 9.1 arcmin. 
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et al. 1989) inside a cooling radius of ~ 1-2 arcsec 
(Sarazin and McNamara 1997). Spectral analysis 
of both the EPIC and RGS data show the, by now 
common, absence of gas cooler than ~ 1/3 of the 
asymptotic cluster temperature of 3.5 keV. How- 
ever we now know that most of the cool gas has a 
different redshift compared to the gas hotter than 
2 keV. Thus for the kT ~ 3.5 keV gas to “cool” and 
become the kT < 2 keV gas, it must also change 
its redshift, which seems rather unlikely. Thus, in 
this cluster, we can reject the idea of cooling from 
our dynamical data. This decomposition of the 
cluster emission argues that most of the mass in 
Abell 2597 resides in the higher redshift, hot com- 
ponent because we see little emissivity from the 
cool compact component, in the RGS data. How- 
ever, a large fraction of all the galaxy reashifts 
are consistent with the lower velocity cool compo- 
nent. It thus seems rather likely that Abell 2597 
has not reached equilibrium between the distribu- 
tion of galaxies and gas. While it is, at present, 
unknown how often “cool cores” have a different 
velocity from the surrounding hot gas, the effect 
we see in our data is rather subtle and may have 
been missed in previous analyses of the RGS clus- 
ter data (e.g. Peterson et al 2003). 

6. Conclusion 

We detect for the first time significant velocity 
differences between different X-ray emitting ions 
in the core of a rich cluster. The lower tempera- 
ture gas is moving at ~ —1320 km s -1 with re- 
spect to the hotter component. The cool compo- 
nents redshift is consistent with the optical and 
radio emitting gas near the central cD and one of 
the two galaxy velocity clumps. The hot compo- 
nents velocity is consistent with the higher redshift 
clump of galaxies. This agreement between the X- 
ray and optical velocities argues that Abell 2597 
is a merger more or less into the plane of the sky. 
The relative emission measures of the cool and hot 
components shows that the bulk of the X-ray emis- 
sivity is associated with the higher velocity com- 
ponent. The different velocities between the cool 
and hot components argues that this object is not 
a cooling flow (in the sense that cooling of the hot 
gas produces the cool gas) and that a more likely 
origin of the cool gas is a merging group. It seems 
less likely, at least in this cluster, that the velocity 
shear is due to the effects of an active galaxy. 


Further analysis of more XMM-Newton RGS 
data and the higher signal to noise and resolution 
data from Astro-E2 will show whether Abell 2597 
is unique. 

Based on observations obtained with XMM- 
Newton , an ESA science mission with instruments 
and contributions directly funded by ESA Mem- 
ber States and the USA (NASA). The Digitized 
Sky Survey was produced at the Space Telescope 
Science Institute under U.S. Government grant 
NAG W-2166. The images of these surveys are 
based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. The plates were 
processed into the present compressed digital form 
with the permission of these institutions. The 
rgsxsrc model was authored by A. Rasmussen. We 
thank J. Peterson for the use of his Monte Carlo 
fitting code. 
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